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Abstract 


The aim of the current paper is to construct the shifted fractional-order 
Jacobi functions (SFJFs) based on the Jacobi polynomials to numerically 
solve the fractional-order pantograph differential equations. To achieve this 
purpose, first the operational matrices of integration, product, and panto- 
graph, related to the fractional-order basis, are derived (operational matrix 
of integration is derived in Riemann—Liouville fractional sense). Then, these 
matrices are utilized to reduce the main problem to a set of algebraic equa- 
tions. Finally, the reliability and efficiency of the proposed scheme are demon- 
strated by some numerical experiments. Also, some theorems are presented 
on existence of solution of the problem under study and convergence of our 
method. 


Keywords: Fractional pantograph differential equation; Fractional-order 
Jacobi functions; Operational matrices; Caputo derivative; Riemann—Liouville 
integral. 


1 Introduction 


In recent decades, fractional calculus and fractional differential equations 
have attracted the interest of many mathematicians and researchers. Frac- 
tional differential equations appear in various fields of sciences and engineer- 
ing such as visco-elastic materials [1], propagation of spherical flames [16], 
electromagnetism [7], fluid mechanics [13], and continuum and statistical 
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mechanics [15]. The fractional nature of this class of equations makes their 
solving very difficult. For this reason, many researchers and authors have 
tried to generalize the existing methods to simply apply and numerically 
or analytically solve the fractional-order functional equations. These meth- 
ods include Spectral method [5,6], Tau method [2], Muntz—Legendre Tau 
Method [17], Collocation method [27], and many others. In this regard, valu- 
able works have been presented, for example ,Dehghan and et al. applied 
the Jacobi polynomials to solve fractional variational problems [4]. In [8], 
the authors utilized new fractional-order Legendre functions to solve weakly 
singular Volterra integral equations. Esmail and et al. in [9] presented the 
collocation Muntz method for the solution of fractional differential equations. 
Kazem generalized the Jacobi integral operational matrix to solve a class of 
nonlinear fractional differential equations [11]. The authors in [12] utilized 
the Tau method together with the d-dimensional orthogonal functions to 
solve partial differential equations. 


In this paper, we consider fractional—-order pantograph differential equa- 
tions as follows: 


l 
D7 u(t) = au(t) + S- bi; DV u(qt) + f(t), m—-1<y<m, te (0,1, 
i=1 
u)(0) = dy, j=0,1,...,m—1, 
(1) 
where a,b; € R, OS % < y < mi = 1,2,...,1,0 < q < 1, D is the 
fractional derivative in the Caputo sense, and u(t) is an unknown function. 


The pantograph equation is one of the kinds of delay differential equa- 
tions and has many applications in electro-dynamic and biology; see [14]. 
Several authors have solved the pantograph differential equations of integer- 
order such as, Jacobi operational method [3], Chebyshev polynomials [21], 
Bernoullie polynomials [23], variational iteration method [26], and so on. 
But there exist few methods applied to numerical solution of pantograph 
differential equations of fractional-order. For example, we can mention the 
Hermite wavelet method [20], spectral-collocation method [24], and Legen- 
dre multiwavelet collocation method [25]. This motivates us to present a low 
cost algorithm for solving this kind of fractional-order equations. To this 
end, we define the shifted fractional-order Jacobi functions based on shifted 
Jacobi polynomial on [0,1]. Then, the operational matrices of fractional in- 
tegration, product, and pantograph are constructed. The resultant matrices 
are utilized to approximate the various terms in equation (1). Finally, the 
main problem is converted to an algebraic equation, which is collocated at 
the roots of PaO). Consequently, an algebraic system is achieved, which 
solving it leads to determine the unknown coefficients vector and thereupon 
an approximate solution is obtained. It must be mentioned that the resultant 
nonlinear algebraic system (corresponding to a nonlinear equation) will be 
solved by the well-known Newton iteration method. 


An operational approach for solving fractional pantograph differential equation 39 


The current paper is organized as follows: Section 1 contains the introduc- 
tion; in Section 2, some basic definitions of fractional calculus are introduced; 
the study of existence and uniqueness of solution of (1) is presented in Section 
3; the shifted fractional-order Jacobi functions (SFJFs) and their properties 
are given in Section 4; in Section 5, the Jacobi operational of fractional 
integration, product, and pantograph are derived; the proposed method is 
described to solve the fractional-order pantograph differential equations in 
Section 6; in Section 7, the convergence of proposed approach is studied and 
an error bound is obtained in the weighted-Jacobi Sobolev space; In section 
8, our numerical results are reported. For this purpose, several examples are 
presented; finally, a conclusion is presented in Section 9. 


2 Basic definitions of fractional calculus 


Some basic definitions and properties of fractional calculus theory are pre- 
sented in this section, which are used in this paper. 


Definition 1. The Riemann—Liouville fractional integral operator of order 
y > 0 is defined as follows: 


t 
DMu(t) = af (t—s)?~*u(s)ds, y>0, t>0, 
0 


Definition 2. The Caputo fractional derivative operator of order y > 0 is 
given by 


1 t q™ 
Dru(t) =I" ID" u(t) = i (¢ —s)™"7-'__u(s)ds, t>0, 
Y) Jo 
where m = [7]. 


For the Riemann—Liouville fractional integral and Caputo fractional deriva- 
tive operators, we have 


1. J Lad) = 1 eu), 

2. I7(Ayur(t) + Ague(t)) = AI? uz (t) + AgI7ue(t), 
T(v + 1) 
I(iy+v+1) 

4, D’ I u(t) = u(8), 


3.1%” = gt’, w>—i, 


’ 


m—1 a (0+) 
i! 


5. YD u(t) = u(t) — 


xr, 
i=0 
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6. Die” = 0, Y> VY, 
ss a ree, otherwise, 


7. DYX=0, 


where Aj, A2, and A are constant. 


3 Existence and uniqueness of solutions of 
fractional-order pantograph differential equations 


In this section, the existence of solutions of (1) is established by using the 
fixed point theorem. 

Let Y be a Banach space and let C(J,Y) be the Banach space of continu- 
ous y(t) with y(t) € Y and t € J = [0,1] and ||y|| = maxze, |y(t)|. Moreover, 
B,(y, Y) represents the closed ball with center at y and radius r in Y. It is 
seen that (1) is equivalent to the following integral equation: 


u(t) =g(t) + cou. (t — s)?~1u(s)ds 
3 Hs y-%-1 ; 
° 2 I(y-%) [ =e) u(qis)ds 


(ae i ae 
+o f (t — 8s)?“ f(s)ds 


where 
m—-1 n m_—1 
dp tk+7 - k 
th= dy thtY-%, m—-1< 4% < mi, 
g(t) d. Tet y+H) ce vi) ja k 4 WS a 
6=1,2,..., 1. 


Theorem 1. [f (ja//[(y+1) +1 bo/Z) < 1/2, where by = maxi<ici {bi} and 
Z = minicicn{l(y-y+)}, then, the fractional-order pantograph equation 
(1) has a unique solution. 


Proof. Let W = C(J,Y). A mapping Vu(t) : W > W is defined by 
t 
a 
Wu(t) =g(t ‘rad t— s)7~1u(s)ds 
(t) =g(t) re : (t — 8)” *u(s) 
t 
al (t — s)?~%~1u(qs)ds 
i) 


+ oT 


l 

i=1 ie 

tr. 7 5a 
aay f, bs oas 
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It should be shown that W has a fixed point and that this fixed point is a 
solution of (1). Set r > 2 (G+ F/I(y + 1)), where G = maxiez ||g(t)|| 
and F = max;e¢,||f(t)||. Then, it can be shown that UB, C B,, where 
B, = {ye W| || y ||< r}. So, one has 


|| Yule) I <I gf) I+ a i Coie ON: 
l t 
Sages fh OI Meee) Ha 
1 


tie [ -9" f(s) || ds 
I 


ae F 
<Gtr 
i, y¥-wl-%) WO) 


la| F 
<G ! ! 
(al ZI Tesh 


<r. 


Thus, © maps B, into itself. Now, for ui (t), u2(t) € W, one has, 


|a| 


t 
|| Yur (t) — Yua(t) Il < | (t—s)7~* || ui(s) — ua(s) || ds 


~ Ty) 
1 
1 [bs | % [7 ae ial u s)—u se E 
Pe ageai, (t— s)7~%™ || ur (gis) — u2(gis) || d 
|a| _l bo 
siya ey eat) 


= lal , 10 ui(s) — ua(s 
(eo +t ZB) e@- ee I. 


As ja|/T'(y +1) +1 bo /Z < 1/2, the mapping VW is a contraction and therefore 
there exists a unique fixed point u(t) € B, such that Vu(t) = u(t). 


4 Fractional-order Jacobi functions 


In this section, first we recall the definitions and some useful properties of 
the orthogonal Jacobi polynomials in the interval [0,1]. Then, we define 
the fractional-order Jacobi functions to obtain the operational matrices of 
fractional integration, product, and pantograph. 
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4.1 Shifted Jacobi polynomials 


The shifted Jacobi polynomials are defined on the interval [0,1], with 
weighted function w'%?)(a#) = (1 — 2)* 2%, and can be determined with 
the following recurrence formula [22]: 


NN al aS, 


2 
E(a ,8,t) P, Pe a), t=1,2,..., . 
where 
,  (+at+fh+1)(a? — B?) 
Alo, 8.1) = se GtatB+ Gitar By’ 
, (+at+B4+2)(%+a+8+1) 
Bla, B, i) 204+DG+e+6+1) 
,_ +a) + D(Q+a+B+2) 
and 
POM) =1, Pa) = SEF Ree 1) 4 oF 


The orthogonality relation of Jacobi polynomials is defined as follows: 
[ ej Pay PPD) wy?) () deo = hy i, 


where 
Tktat+1l(k+6+4+1) 


Qk+a+8+ 1k! T(k+at+6+1) 


The analytic form of Jacobi polynomials is given by 


hy = 


POP) (2) a5. (1 "T+ B+ DIG tk +046 412" 


24 Tk+ B+ G+at+B+1) GWM ” 


A continuous function u(t), in the interval [0,1], can be expressed in terms 
of shifted Jacobi polynomials as follows: 


u(z) = > C; PI~)(a), (3) 


where the coefficients C’; are given by 


1 1 
ae ~f u(t) P(e) w(x) de, 7 =0,1,--. . 
Jj JO 
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In practice, only the first (N + 1)-terms of shifted Jacobi polynomials are 
considered. Therefore, one has 


N 
un(x) = >_C; Pi) (2) = OT (2) C=CT O(a), (4) 
j=0 


where the vectors C and ®(t) are given by 


C=[Co,C1,...,Cn]?, B(x) = [PI (ew), PI (a), ..., PR (ae). 


(5) 
Also, the following auxiliary relation is obtained from (2). 
(2,8) ! (2,8) A(a, 8,t) p(a,8) 
Qe P = P' _ P' 
oe ) Bia, B,i) +} (x) B(a,8,i) * (x) 
+ Fle. B1) p(o.8) (5) i;-1,2 N 
A a eee . 
Bla, B,i) * oe ea 


Two other properties of the shifted Jacobi polynomials are defined below: 


P= (ay (*F%), 


(7) 


EPP (a) Pn titatBtl) parses 
dat T(nt+a+B+1) 


(x). 
4.2 Shifted fractional-order Jacobi functions 


We define the shifted fractional-order Jacobi functions (SFJFs) by trans- 
formation x to t? (a > 0) based on Jacobi polynomials (2) in the inter- 


val [0,1]. For convenience, the fractional-order Jacobi functions, pie (te), 
i = 0,1,2,..., are denoted by ploPie) (t). It is clear, the classic Jacobi 


polynomials are resulted for o = 1. The weight function is w'%*7)(t) = 
ot7(8+1)-1(] — t7)%, Moreover, the relations (2)-(7) are converted to: 


PY? (1) =Ala, Bd) PIP (t) + (t? — 1) Bla, B, 4) PI (8) 


(8) 
— E(a, B,i) PP”) (2), a ae 


0 = 
PioB.2) (t) bes i; Pk.) (4) = a+ s + a B 


where coefficients A(a, 6,7), B(a, 8,1), and E(a, 6,2) are same in (2). 


1 
| PLP) POP (t) wlOP)(t) dt = he jx, 
0 
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(-1)' "TG + 6+19TG+k+a+8+1)t* 
TkK+6+1Gta+64+1) G—k)! kl’ 


(9) 


where 


a @,8,0 
ee ~{ u(t) PP) (4) we) (2) dt, 7 =0,1,---, 
Jj JO 


N 
un(t) = S10; PLP (t) = BIT) C = CT O(N), 
j=0 


where the vectors C' and ®(7)(t) are given by 


C=[05,CyigOy, BOG =[POP I ®, POPP On, POO el. 


(10) 
o pla,B,c) 1 (a,8,c) A(a, B,1) (a,B,c) 
27 P. (=. 7 poe e ' 
E(a, B,1) (a,8,c) 
+ Basa Perse c(t) tH 12.2. GN, 
gO 


5 Fractional-order Jacobi operational matrices of 
integration, product, and pantograph 


The objective of this section is to derive the operational matrices of fractional 
integration, product, and pantograph related to the SFJFs. The fractional- 
order operational matrices of product and pantograph are utilized for ap- 
proximating the product of the vectors ®(7)(t) and ®(7)7(t) and the terms, 
which include the delay in equations under study, respectively. First, some 
properties of the shifted fractional-order Jacobi functions are presented as 
follows. 


Lemma 1 (see [19]). The ith shifted fractional-order Jacobi function, p87) (t), 
t € [0,1], can be obtained in the form 


n 


Pern 
k=0 
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where coefficients a are given as 


af “ae Gace Gar 


Lemma 2. [fn EN and pn, then we have 


1 

/ tP P&P) (4) ay (8) (4) dt 
0 

_ > (-1)"-*§ T(nt+ 641) T(n4 

<p OTR + B+) Tmt+a+B 


k+a+8+1)T(k+8+p/o+1) Pat) 
1) T(kK+a+6+p/o+2) (n—k)! k! 


Proof. The lemma can be easily proved by multiplying relation (9) in t? and 
then integrating of the resulting equation. 


Now, the Jacobi operational matrix of fractional integration is generally 
derived. 


Theorem 2. Suppose that ®7) (t) is the shifted fractional-order Jacobi vector 
in equation (10). Then, we have 


Ta) (t) ~ P™ 6 (t), 


where I‘) is the Riemann-Liouville fractional integral operator of order v 
and P\) is the (N +1) x (N +1) operational matrix of fractional integration 
which defined by 


(0,0) (0,1) m(0, N) 
om m(1,0) (1,1) m(1,.N) | 
n(N,0) (N,1) ... (N,N) 
where ; 
=> Wey 2S 01 7 1 (12) 
k=0 


Ae 
= 
Q 

+ 
& 
pe 
a 
+ 
= 


!(@—DW(ol+v +1) 


Proof. By using the analytical form (9) and the properties of operator I”, we 
have 


tat B+1)P(ol+1) +7 
1) P(ol+v+41) i! (i—D! 


(13) 


a 2 (-1)'-! PG + BED) TG 
1” ph BD (4) dt = > ( ai , . e aI Ut (2 
1=0 T T 
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Now, ¢?/+” can be approximated in terms of the SFJFs as follows: 


ett a Ss P PloF.2)( (t), 


where 


1 f sie 
ang = 5 | gee PONG) ah Piel (a) dt, 


According to Lemma 2, one has 


al 
| polity re a(S) we) (£) dt 
0 


_ACDIRH* Mat) TG4 
2 Te+e+DTG+a 


B 


co 
9 |+ 
DR 
is 
ne, 
q 
XL 
= 
oS 
& 


Bo D 
Therefore, equation (13) is written as follows: 


pe plate? (t dt 


N a . 
re (1) TG + B41 TG+64+D)TG+l+a 
yp» hj TG tot+6+1)TU+84+1)Téta+B 


jo \ kao t +a )Tol+v+1)UuG-D! 
1. (-1-* TU th ta+Bt1) Pl+k+B+v/o4+1)) p(a8,0) 
x. Taree TC k+a+B+v/o+ G nara ) 
N 
=> ali) POPE) 
j=0 


(14) 
where 7(i,7) are given by equation (12). Accordingly, equation (14) can be 
rewritten as the following vector form: 


1” PP?) (t) ~ [x (i, 0), w(i, 1), ---, ti, N)] POW), §=0,1,...,N. 


This leads to the desired result. 


The following lemmas are useful to obtain the fractional operational ma- 
trix of product. 


Lemma 3 (sce [{19]). Jf Perey and Peer”) (t) are respectively jth and 
kth SFJFs, then the product of PONG) and pla) (t) can be written as 


Q's B,o) sk) yor 
Qik -S a t 


where coefficients A; @*) are determined as follows: 
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Ifj>k: 
r=0,1,...,79 +4, 
if r>J, then 
ik ; k 
Ane wis hae 1 i 
else 
ry = min{r, k}, 
ik Tr j k 
A) = ee Hs 
end. 
Ifj<k: 
r=0,1,...,7 +4, 
if r<g, then 
ry = min{r, jf, 
j,k r } k 
ay v= a8 ae vf ' 
else 
ro = min{r, k}, 
j,k r j k 
xe i ay ge vf : 
end. 


The quantities y™ and oy, are introduced, respectively, for poe 7) (t) and 


pee (t) based on the Lemma 1. 


Lemma 4. [f pA) (t), pee (t), and plmae) (t) are respectively i-, j-, 
and kth SFJFs, then, 


TF  (1y-e) TE et 


= ) 
a9 T(l+8+1) Tita 


where ae are computed by Lemma 38. 


Now, a general formula is presented for finding the (N + 1) x (N +1) 
operational matrix of fractional product U whenever, 


8 (4) FT (4) U~U & (2), (15) 


and U is a (N + 1) given vector. 


Theorem 3 (see [3]). The entries of the matrix U in (15) are computed as 
follows: 


N 

2 1 ; 

Ui = y Ui dijk, j,k =0,1,...,N, 
1=0 
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where qijr are introduced by Lemma 4 and U; are the components of the vector 
U in (15). 

Corollary 1. If u(t) ~ UT & (t) = 67 (4) U and &) (t) BOT (4) Ux 
U 6(7)(t), where U is an (N +1) vector and U is the (N +1) x (N +1) 
operational matrix of fractional product, corresponding to the vector U, it is 


easily shown that, - 
u(t) ~ UT 0 Ot). 


To approximate the terms including the delay, the fractional-order Jacobi 
operational matrix of pantograph is derived. For this end, we apply the 
auxiliary (11). 

Replacing t with gt in definition (10), leads to 


BN (gt) = [POP (ge), BPP Gt) ng Pee a] 


The following theorem expresses the relation between vectors ®(7)(qt) and 
#7) (t). 
Theorem 4. The relation between the delay vector ®(~) (qt), 0 <q <1, and 
vector ®(7)(t) can be expressed as follows: 

Ho) (tq) = L(qc) #7) (t), 
where Lig,c) is an (N +1) x (N +1) lower triangular matrix called fractional 
operational matrix of pantograph. 


Proof. Replacing ¢t with gt in (8) and rearranging leads to 


PEP (qt) =A(a, Bi) PP (gt) + (2q7t? — 1) Bla, 8,4) PLP (at) 


i4l 
— E(a, B,i) PCP (qt), 1 =1,2,... - 
(16) 
The ith delay SFJF , PioP-) (ge), can be expressed in terms of the SFJFs 
themselves as follows: 


APNG) = The BPM, EHO AD 
k=0 
From the first few polynomials in (17), it easily obtains 


loo=1, logs =0, ¢=1,2,...,N, Lo =(8+1)@’-1), hi=d’. 
(18) 
It is clearly found that, 
ag, 220... (19) 


Using auxiliary (11) and substituting (17) into equation (16), we have 
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i+1 


ying eee |= 
k=0 
ee, ge ue , A(a, Bk), (a,8,0) pie eer 
= > 4 { (Alsi +4 Bla, B,i)(1— Bla, Bk) Pot (t) - B(a,8,1)) P. (t) 
: 1 ; E(a, 8, k) (a,B,o) 
+ 4° B(a, B,1) Conn ' Ble, Bk) P3 )} 


ft 
— E(a, B, i) » li-ayk es 
k=0 
(20) 


Equating the coefficients of Pohl) on both sides of equation (20), for 
k =0,1,...,7%, leads to 


biti,0 =| Ala B, 1) 


; . wk 2at1)(6+1) 
+4q B(a,8,i)(1 ot Gee) 


- Blo, .i)} lio — E(a, 8,2) li-1,0, i= 12 eieg dN 1, 


ae ={ Ala, 8.9) +47 Bla, Bi) FES - Bla + B+ »)} . 


: E(a,8,k +1) : | 
cra a ka + Bq, BI) kins} q° B(a, 6, i) 


— Cla, Bt) li-tk,, k=1,2,...,4-1, 


A } ; 
—e (4c, Bi) + q°B(a,B, (1 — = ~ Ba, B, i) fd 


PO in, PHB A, 
Bia, B,i— 1) 


The starting values are given by equations, (18)—(19). Finally, for the pro- 
portional delay g, we have 


G7) (gt) = [Pk-?”) (gt), Ph? (tq), ..., PP (tg) |? 


a,8,0) 

1 0 0 ... 0 0 Fo : t 
he OF © «as 0 6 |) |PeeOG 
Po. 0 DO | | PROm 


loo loa 


In—1,0 lw—11 lw—12--. GX? 0 P ofie) (4) 


N 
l l l Yecal ip Qe a,B,o 
N,O Nl N,2 N,N-1 4 pK B (t) 
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= Liao) 8°) (a). 


Corollary 2. The operational matrices of fractional integration and product 
are converted to the following forms for the delay q. 


P’@O (gt) ~ Liao) P™ 8 (8), 


B) (gt) BOT (gt) Cy Lig) FOO), F= beac, 


where the matrices P) and L(q,0) were defined by Theorems 2 and 4 and 


F is the fractional operational matrix of product corresponding to the vector 
F= Lig a C. 


6 Description of method 


Consider the problem given in (1) again. We expand D7u(t) by means of the 
SFJF's as follows: 
Dru(t) ~ BT (£)C. (21) 


According the properties of Caputo fractional derivative, we have 
m—-1 d 
~ IPT ok sk 
u(t) ~ POT (t)C + > ae 
=0 
~ OT (A)POTC+ SOOT EF, m-l<y<m, 


(22) 


where F is a known vector and 3. «tk is approximated as ®()" F. To 


approximate the Dviu(t),0 < % <7 < m,i = 1,2,...,1, the following 
procedure can be pursued. 


D u(t) = D!~-% D¥ u(t) ~ BT (4)C. (23) 


Using the Caputo derivative, we get 


m—1 
d 
Diiult) ~ IY GET (4 k k—i 
u(t wet SS paesicay 


~ GOT) PO-WTICYOOTR, OS y%<7<m, i=1,2,...,1 
(24) 
Now, by applying equations (21)—(24) to the terms containing the delay, we 
get 


An operational approach for solving fractional pantograph differential equation 51 


u(qt) > GOT H)LE oP OTC 7 oT (LE ae 


(a, 
D¥u(gt) ~ BOT OLE PO wWTC+ EO MLE A, $= 1,2,...51, 
2 ~ TTT AHl(o)T —_ p(y)T w pHlo)T T 
u(t) ~U* U8" (t), U=PV°C+F, Du(qt)2 © via 


where U is the fractional operational matrix of product corresponding to the 
vector U. By substituting the expressions (21)—(25) into (1), the following 
algebraic equations is resulted. 


BOT (HC —aBOT()(POTC + F) — S°b,67 (t)(LE, PO WTC 
i=l 


+ Lé,,0) Fs) — f(t) © 0. 
In this way, the main equation is converted into an algebraic equation. The 
resultant algebraic equation is collocated at (N +1) roots of the (N + 1)th 
shifted Jacobi polynomials on the interval (0,1). By solving the system of 
generated algebraic equations, the unknown vector C’ can be determined. 
Thus, an approximate solution can be acquired for various values of a and £ 
parameters using equation (22). 


Corollary 3. The nonlinear terms in the problems under study are approx- 
imated as a matrix form. For this reason, it is necessary that the initial 
conditions of the problems to be approximated in terms of fractional-order 
Jacobi basis. 


7 Convergence analysis and error estimation 


In this section, the upper bounds are presented for the errors of the residual 
function, proposed algorithm, and operational matrix of fractional integra- 
tion, P‘). It is seen that the upper bounds tend to zero when the number of 
the SFJFs increases. First the following definitions and theorems are stated 
for 0 = 1. 

The set of all algebraic polynomials of degree at most N is denoted by Py. 
The orthogonal projection Py,«,gL? (.,s)(J) > Pn, J = [0,1], is considered 
for u(t) € L? ...)(J) and defined by 


(Pn,a,ev —v, >) =9, for all 6 € Py. 
The Jacobi-weighted space is introduced as 
Fita.py(J) = {v| v is measurable and ||v||,u(a0) < co}, reéN, 


equipped with the following norm and seminorm, 
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r z 
ibson (> | Fv [osreen ) a eee eee 


k=0 


Theorem 5 (see [10]). For any v € Fr .a.3)(J), rE NandO <u <r, we 
have : 
lv — Py,a,6Ull wie) < C(N(N + at B))F lol, (o.8)» 


where c is a positive constant independent of N, a, and B. 


Lemma 5. [fn EN, v€ Fr.a.4)(J), and D" v(t) is the derivative of v(t) of 
order n, then 


ur 


2 
||[D°v — D®Py,0,6Ull pws) < c (w +n)(N+n+a+ 8)) |v |, ewan) « 


Proof. By derivation of the shifted Jacobi polynomials and applying Theorem 
5, the desired result is acquired. 


Corollary 4. Forv € Jiiaa)(J) and0 <q <1, we have 


Ilv(qt) — Prv,a,82(4t) Ilene) < e(QN(N +a+ B))°= |v(gt)|nwe- (26) 


Theorem 6. If y > 0, y € R, and D v(t) denotes the Caputo fractional 
derivative of u(t) © Tiva,a,(J) of order y, then 


por 
2 


rol 


Dv — D’Prve as) < 
|DPv — DY Prv,a,6%lluwe < Tana 4) 


(w m)(N +m ++ 8)) 
x |e | ayler) m= [yl]. 


Proof. Using properties of the Riemann-—Liouville fractional integral operator 
leads to 


D v(t) =I" ID" v(t) = aes) 7 (¢ — s)™—7—1 D™v(s)ds 
1 m—y—-1 x D™y 
=i D™v(t) 


9 


where the star symbol denotes the convolution of two functions t”~7~! and 
D™ v(t). So, by means of Lemma 5, we have 


||D%v — D? Py, 0,8 lyw(oa = |Z" 7(D™0 — D™ Py, 0,60) Ilu,w6a®) 
¢m—-y-1 
Toa) * (D™u — D™ Py a.pv) 


p,w(o,8) 


1 
< 


pb™ D™ - a 
Gn Fm 29) IBY — DOP r..8Y)Ilawto 
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bor 
/ 


Cc 2 
<= | (N N NT ails 
ren (WW HMO +m a+8)) © Llp wian 


Now, we consider equation (1). If uy(t) is the approximate solution to 
u(t), then the residual function is as follows: 


l 
Hy (t) = D7un(t) — aun (t) — 5b; DY un (ait) — f(t). (27) 


i=1 
The next lemma estimates a bound for the residual function (27). 


Lemma 6. Consider the residual function given in (27). A bound for it can 
be estimated as follows: 


wor 
t 
Cc 2 


T(m—y+1) 


ta — 1 («C20 + my (a +m +a+3)) 


[F(a <O*( (w+ mj(v+m a+ 8)) 


wor 
2 


M- 


i=1 


+ lale(W(N Fat p's), 


aX trey |u*) (t) |, ay(a.e), =0,1,...,m, mj—1<%< mj, and 


Proof. By using Theorems 5 and 6 and (26), the above inequality can be 
resulted. 


The following theorem presents an error bound for the proposed method. 


Theorem 7. If uy(t) is the Jacobi approximate solution to u(t) and ex (t) = 
u(t) — un(t) ts the error function, then the error bound can be estimated as 
follows: 


ur 
2 


(1/Ja| + 1)c 


lex (4) II pw.) (Tmt i) (Cr + my +m a+ 8) 


ur 
2 


+ yf see = eee (aw +m;)(N+m,+at 9)) 


oer 


Proof. By subtracting equation (26) from equation (1), we get 
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I 
Dien(t) = aen(t) + _ b; DV ey (git) — Hy/(t). 


i=l 


Now, applying Theorems 5 and 6, equation (26), and Lemma 6 leads to the 
desired result. 


Now, we consider the operational matrix of fractional integration P™ in 
Theorem 2. We define the error vector FE, as follows: 


E, = 18) (t) — P 8) (4) = [Eoy, Fiy,..-, Env)", 


where 


N 
Bp =P PPO) — SPO, POP" ®), k= 0,1,...,N. 
j=0 
The following theorem presents an upper bound of Riemann—Liouville frac- 
tional integral operator. 


Theorem 8. If Ex, = Peery. Pe per (t), then an error 
bound of Riemann-Liouville fractional integral operator of order v can be 
expressed by 


Ex vl, yas) SO(N(N +o + 8))B-" 


kk OT (oj + WP (oi + DE(oj tv + 20 (oi +v +2) 
x im [p+ ilote y+ it) +y— 9. al (aero) 


j=0 i=0 


x Bi+jtr+B+2(v Noriarreyh, 


where ~y$h) were introduced by Lemma 1 and B(u,v) is the well-known Beta 


j 
function. 


Proof. By using Lemma 1 and the property of Riemann—Liouville fractional 
integral, we have 
: T(oj +1) 
POPP ye. Sy gee, 
& (1) du T(oj+v+1) 
By setting Z = 1” P\?")(t), we get 


2 = 2 
IZ) (ae) = WOE ZI catretr) 


_ | V(ojg + DI (of + +2) 
T(oj tv+1)(oj +v—1r +2) 


2 


tri ter 


pywlotr.B+r) 


oy MY D(aj +E (oi+ DE (Oj + + 2)0 (ci + + 2) 
im ty+ 1)(oitv+ (oj t+v—r+2)[(oit+v—r +2) 
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aL 
x | to tj+r+B8+1)+2(v—-r) la = reyetratl 
0 


kk yD (oj + IL (oi + IL (oj +v + 20 (oi + +2) 
= » lie 


a 2) jtv+Di(oitv+Di(oj+u—r+2(citv—r+2) 


L 
x | gititrtB+2v—r)/o cy = neta, 
0 


Therefore, the theorem is proved. 


8 Numerical examples 


In this section, the method presented in this paper is applied to solve some 
examples taken from [18] to demonstrate the efficiency and reliability of the 
proposed algorithm. As a result, a comparison is presented. 


Example 1. First, consider the following fractional pantograph differen- 
tial equation from [10], 


3 t 
Du = que) +u(5) — +2, 1<7<2, (28) 


with conditions, 
u(0)=0, ~u'(0) =0. 


The exact solution is u(t) = ¢? in the case y = 2. We first approximate the 
unknown function and its derivative of fractional order y as follows: 


t 
ur(t) ~ BOF (4) C, u(t) ~ BOT (t) PMY C, u(5) » )T (t) Liao) POT C, 


Substituting the approximations into (28) leads to the following algebraic 
equation: 
3 


FT raha POT C_ OT (t) Li o) POT C+1?-2 0. (29) 


By choosing N = 4, equation (29) is collocated at roots of pio) (t). Then, by 
solving the resultant algebraic system, the unknown vector C' can be deter- 
mined for some values of a and 3 parameters. Table 1 displays the absolute 
errors at some arbitrary points for a = 8 = 0, y = 2, 0 = 0.5,1,1.5,2, and 
N =4. Table 1 shows that the results of our algorithm are more accurate 
than those which reported by [18]. In Figure 1, the approximate solutions are 
compared to the exact solution for a= 6 =0, N = 4, y = 1.3,1.5, 1.7, 1.9, 2, 
and 0 = 1. From these results, it is seen that the approximate solutions 
converge to the exact solution when y approaches 2. 
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Table 1: Errors for a = 8 = 0, N =4, and various values of o for Example 1 


ti c=0.5 c=1 c=15 c=2 oc =0.5 c=1 o=15 
in [18] in [18] in [18] 
0.0 5.3967 x 1072 1.0616 x 1075 1.5319 x 1073 6.4935 x 1075 5.02 x 1072 1.64 x 1071 1.33 x 107-1 
0.2 | 9.5185 x 1074 5.6130 x 1075 4.2342 x 1074 1.1404 x 1074 = — — 
0.4 1.8965 x 1073 9.4037 x 1074 2.0913 x 1073 9.7976 x 10~4 3.08 x 1072 6.51 x 1072 1.47% 107° 
0.6 4.4029 x 1073 4.7621 x 1073 3.4555 x 1073 4.8110 x 1073 == —— — 
0.8 | 1.4375 x 1072 1.5160 x 107? 1.5934 x 107? 1.5206 x 1072 1.59 x 1071 2.60 x 1072 3.78 x 1071 
1.0 3.8299 x 1072 3.7313 x 107? 3.1948 x 1072 1.5206 x 1072 2c ae — 


T T T T T T 
0 0.2 0.4 0.6 0.8 1 


— exact solution —~y=1.9 °""" y=al.7 —-~ y=15 —— y=13 


Figure 1: 
Example 1. 


Plots of u(t) for N = 4, 0 = 1, y = 1.3,1.5,1.7,1.9 and exact solution for 


Example 2. Consider the following linear fractional-order pantograph dif- 
ferential equation of [18], 


Dvu(t) = —u(t) +0.lu(St) 40.5D"u( st) + (0.32405) 0 Rien), 


where 0 < y < 1. With condition u(0) = 0 and the exact solution is u(t) = 
texp(—t). By using the fundamental matrices obtained in Section 5, the 
following algebraic equation is resulted: 


BOT (t) C+ HOT) POT C— 0.107 (t) Lira 5)” POT C 
4, 


~ 0.587 (t) Lia)’ C — (0.32t — 0.5) exp(—0.8t) — exp(—t) ~ 0. 
(30) 
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By choosing N = 10, equation (30) is collocated at roots of Pho) (x). Then, 
by solving the resultant algebraic system, the unknown vector C’ can be 
determined for some values of a and ( parameters. Table 2 displays the 
maximum absolute errors for various values of a and § parameters, N = 10, 
and y =a =1. The data of Table 2 shows that the approximate and exact 
solutions are in a good agreement. In Table 3, the values of absolute error are 
compared to the values reported in [18] at some arbitrary points for N = 10, 
a = B = 1/3, and y = o = 1. It is seen that the errors of our method are 
lesser than the errors of methods applied in [18] on the interval [0,1]. Also, 
the absolute errors are listed in Table 4 for N = 10, a= 8 = 1/3, y =o and 
various values of 7 in some arbitrary points within interval [0,1]. A graphical 
comparison between the exact and approximate solutions is depicted in part 
(a) of Figure 2 for values of a = 8 = 1/3, y =o = 1, and N = 10. From 
part (b), it can be shown that the numerical results are well consistent with 
the analytical solution. In addition, the plots of u(t) are shown in part (0) of 
Figure 2 for o = 1 and various values of y. It is seen that the approximate 
solutions converge to the exact solution as y approaches 1. 


Table 2: Maximum absolute errors for various values of a and 8, y = ¢ = 1, and N = 10 


for Example 2 


5.5014 x 1078 
1.5979 x 10712 


— 


—2) | 8.6704 x 10-"4 
3) | 1.1469 x 107” 
2 


(a, 8) Error Abs a, 3) Error Abs 
(0,0) | 2.2652 x 10-7 3,3) | 4.5600 x 10-8 
(1,1) 7.6421 x 10-8 353) 1.1765 x 10-8 
3 
2 


aay 
| les 
NY wits 
Soa eee 


hicowIhy ps “CD>| O 


— | 


Table 3: Comparison of absolute errors with other methods for Example 2 


tj present method Bernoulli Runge-Kutta Variational iteration @-method 
method in [18] | method in [18] method in [18] in [18] 
0.1 | 6.4471 x 10-14 1.98 x 10— 8.68 x 10-7 1.30 x 10~3 4.65 x 10-3 
0.3 4.4701 x 10714 7.78 x 1079 1.90 x 1073 2.63 x 1073 2.57 x 1072 
0.5 | 4.7819 x 10716 6.34 x 1075 2.28 x 1073 2.83 x 1073 4.43 x 1072 
0.7 | 4.5170 x 10714 4.36 x 1075 2.28 x 1073 2.39 x 1073 5.37 x 1072 
0.9 | 6.4785 x 10714 2.80 x 1075 2.03 x 1073 1.64 x 1073 5.35 x 1072 


Example 3. In current example, we consider the following fractional pan- 
tograph differential equation of [18] : 


t t. 1 t 
} D™ 
3) us 


D u(t) = sult) + u( 
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Table 4: Absolute errors for various values of y = 0, a = 8 = = and N = 10 for Example 
2 
tj 7 =0.6 y=0.7 7 =08 7 =0.9 
0.2 | 2.5022 x 10-5 | 1.1429 x 10-4 | 2.6347 x 10-4 | 3.3616 x 10-4 
0.4 | 1.0801 x 10~° | 5.3195 x 10~° | 1.3757 x 10-4 | 1.9924 x 1074 
0.6 | 6.3865 x 10~® | 3.1052 x 107° | 7.9303 x 107° | 1.1741 x 10-4 
0.8 | 4.2784 x 107-© | 1.8443 x 107° | 4.4219 x 107° | 6.5747 x 1075 
1.0 | 5.0278 x 10-© | 2.0332 x 10-® | 4.6987 x 107° | 6.8103 x 107° 


T T T T T T T T T T T 
0 02 04 06 08 1 0 0,002 0.004 (0,006 0.008 0.010 
x x 


— exact solution fory=o=1 © approximate solution — exact solution ~~ y=0.9 °-""* y=08 =" y=07 —~ y=0.6 


Figure 2: Plots of (a) exact and approximate solutions, (b) u(t) for a = B = 5 o=1, 
7 = 0.6, 0.7, 0.8,0.9, and N = 10 for Example 2. 


where 0 < y1 < y < 2, with conditions u(0) = u/(0) = 0. The exact solution 
is u(t) = ¢? for case y = 2 and y, = 1. Substituting the approximations 
introduced in the previous sections into (31), leads to the following algebraic 
equation: 
eT (t) C- SH") POT C_- OT (Et) Lit.) PTC 
[iia aa 3 a 5B (t) Ling) C+? +t-180. 
(32) 
By choosing N = 10, equation (32) is collocated at roots of P&P) (x). Then, 
by solving the resultant algebraic system, the unknown vector C’ can be 
determined for some values of a and ( parameters. Table 5 displays the 
maximum absolute errors for various values of a and @ parameters, y = 2, 
1 = 1,0 =1, and N = 10. The approximate solutions are in good agreement 
with exact solution. Equation (31) was solved in [18] by Bernoulli collocation 
method. In Table 6, the absolute errors are compared to the values reported 


— plat (t) Li e) 
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by [18] at some arbitrary points for N = 10,a=6=2,y7y=2,mn =1, 
o =1. Table 7 shows the maximum absolute errors for a = 6 = 0, y = 1.9, 2, 
1 = 0.9,1, and o = 0.5,1,1.5,2. For the case y = 1.9, 7; = 0.9, and the case 
y = 2, 71 = 1, the best results obtain for o = 0.5 and o = 1, respectively. 
The errors given in [18] are greater than the errors of our method at the 
selected points on the interval [0,0.5]. A graphical comparison between the 
exact and approximate solutions for values of a = 6 = 0, 7, =o = 1, 
y = 1.2,1.4,1.6,1.8, and N = 10 are shown in Figure 3. 


Table 5: Maximum absolute errors for various values of a and 8, y = 2, y1 = 0 = 1, and 
N = 10 for Example 3 


(a, 8) Error Abs (a, 8) Error Abs 

(0, 0) 7.2345 x 10-2 | (3,5) | 1.9150 x 10- 

(1,1) | 1.8407 x 10770 as —#) | 8.7002 x 10~?? 
(—§,—q) | 2.2935 x 107° (2,2) 12207 x10 


Table 6: Comparison of absolute errors with other methods for Example 3 


tj present method Bernoulli Runge-Kutta Variational iteration 0-method 
method in [18] | method in [18] method in [18] in [18] 
0.1 0.0000 3.47 x 10~ 16 1.00 x 1073 1.67 x 10-4 6.10 x 1073 
0.2 1.0000 x 10721 5.55 x 10716 2.02 x 1073 7.15 x 1074 2.58 x 1072 
0.3 | 2.0000 x 10721 3.05 x 10~ 16 3.07 x 10-3 1.73 x 1073 6.47 x 1072 
0.4 0.0000 1.14 x 10715 “17 « 16-3 3.30 x 1073 1.37 x 10-2 
0.5 1.0000 x 10729 1.39 x 10718 5.34 x 1073 5.50 x 1073 2.81 x 1072 


Table 7: Comparison of maximum absolute errors for different values of y, 71, and o for 
Example 3 


ag =0.5 g=1 g=1.5 o=2 
2.2058 x 1077 | 1.0822 x 10~? | 1.0081 x 107? | 1.5358 x 1071 
y=2 n=l 
a =0.5 ag=1 a0=1.5 o=2 
3.4229 x 10-78 | 7.2345 x 10-™ | 6.8343 x 10-9 | 1.1564 x 1077 


Example 4. Consider the following pantograph fractional differential equa- 
tion [18] : 


t t 
5 + 9u( 


D7u(t) = (t) + 4u( y+P-1, O<7<1. 


oy 
6 
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— exact solution —~- y=1.8 ----° y=1.6 —— y=1.2 —:- y=14 


Figure 3: Comparison of exact and approximate solutions for y1 = o = 1, y = 
1.2, 1.4, 1.6, 1.8, and N = 10 for Example 3. 


The initial condition is u(0) = 1. The exact solution is u(t) = 1+ Sa + 
He 7? + Lite? in the case y = 1. By applying the fundamental matrices 


presented in the previous sections, we get the following algebraic equation: 


BOT (tL) C+ {Hr (t BOW GOO) Fp aoe G) Te. 3 POC 
a 
elt (e) LaF} 9{ B17 (4) Lit.) POT C 
+ BT (¢) Lis.) F}-? +120. 
(33) 


where the matrices L, 1) and L, 1g) are pantograph operational matrices 


corresponding to the delays q = 3, 4 and u(0) is approximated as H()T (t) F. 
By choosing N = 10, equation (33) is collocated at roots of PLP) (g), By 
solving the resultant algebraic systems, the unknown vector C' can be deter- 
mined for some values of a and @ parameters. In Table 8, the absolute errors 
of the approximate solutions are listed at some arbitrary points for N = 10, 
a= 8 =7=1, and various values of a. As it is seen, the best results are re- 
lated to the values 0 = 0.5,1. A graphical comparison between the exact and 
approximate solutions is found in part (a) of Figure 4 fora = B =y=o=1. 
The plots of u(t) have been compared in part (b) for y = 0 = 0.75, 0.85, 0.95. 
Also, the error functions of these solutions are depicted in part (c). It is seen 
that the approximate solutions converge to the exact solution as y and o 
approach 1. 
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Table 8: Absolute errors for various values of ¢, a = 8 = y = 1, and N = 10 for Example 


4 
tj o = 0.5 o = 0.75 g=1 o0 =1.5 
0.2 | 1.5100 x 10-™ | 1.1748 x 10-3 | 2.0000 x 10-19 | 4.5227 x 10-1! 
0.4 | 3.5176 x 10715 | 2.3520 x 1073 | 9.0000 x 107!9 | 1.0366 x 107° 
0.6 | 6.5190 x 10715 | 4.5162 x 1073 0.0000 1.9337 x 10~° 
0.8 | 1.0666 x 10714 | 7.3187 x 1073 | 5.0000 x 107!8 | 3.1543 x 107° 
1.0 | 1.6149 x 1074 | 9.6215 x 107 | 3.0000 x 10-18 | 4.7377 x 107° 


(a) 


T T T T T T 
0 0.2 04 0.6 08 1 
x 


—— exact solution © approximate solution 


(c) 


~ 
12+ 
i 
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We vact ” “approximate! y=18 = y=L7 7 Y=L6 


Figure 4: (a) Exact and approximate solutions for a = 8 = y = 0 = 1, and N = 10, 
(b) comparison of u(t) for y = o = 0.75,0.85,0.95,1, (c) maximum absolute errors for 
7 =o = 0.75, 0.85,0.95 for Example 4. 
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Table 9: Absolute errors for various values of o, and N = 10 for Example 5 


t; o = 0.5 o=1 o =1.5 o=2 

0.2 | 2.2499 x 10-38 0.0000 5.0015 x 10-° | 2.4119 x 10-2 
0.4 | 1.3489 x 10747 | 1.0000 x 1072! | 6.8587 x 107° | 7.0770 x 1072 
0.6 | 3.9130 x 10717 0.0000 5.8770 x 107° | 1.6689 x 107 
0.8 | 8.4580 x 10747 | 4.0000 x 10728 | 2.1854 x 1074 | 3.1842 x 107! 
1.0 | 1.5572 x 10726 | 1.0000 x 10729 | 2.3271 x 1074 | 5.5809 x 1072 


Example 5. In this example, we consider the following fractional pantograph 
differential equation of [18] : 


t t 1 t 
Dru(t) = u(t) + D*u(5) + DY u(s) + sD%u(5) — tt — 5t 


where 0 < 71 < y2 < Y < 3. The initial conditions are u(0) = u’(0) = 
u"(0) = 0 and the exact solution is u(t) = t* in case 7, = 1, y2 = 2, and 
y = 3. By using the approximations introduced in the previous sections, we 
get the following algebraic equation: 


ET (4) C- OT (4) PMT C— OT (t) Lis.) PO-wWT CG 
1 


= oo) (t) Lio) po-2)T C- so Les) C (34) 
1 4 
i? + 2? + 2? — 21t 0. 
+ 5 + 3 


By choosing N = 10, equation (34) is collocated at roots of Pho) (x). Then, 
by solving the resultant algebraic system, the unknown vector C is deter- 
mined for some values of a and $ parameters. Table 9 displays the absolute 
errors for a = 8 = 0, y = 3, ye = 2, 1 = 1, and various values of a. In 
Table 10, the values of the approximate solutions are shown at some arbitrary 
points for N = 10,a= 8 =0, y = 3, y2 = 2, 71 = 1, and o = 1. A graphical 
comparison between the exact and approximate solutions is found in part (a) 
of Figure 5 fora = B =0,0 = 7% = 1, y2 = 2, and y = 3. The plots of 
u(t) have been compared in part (b) for y = 2.2,2.4,2.6,2.8,0 = y, = 1, 
and y2 = 2. It is seen that the approximate solutions converge to the exact 
solution as y approaches 3. 
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— exact solution —- y=2.8 —— y=2.2 -"~ y=2.4 —~y=2.6 


63 


Figure 5: Comparison of (a) exact and approximate solutions, (b) u(t) for 
1, and o = 1, and y = 2.2, 2.4, 2.6, 2.8 


a=(6=0, and N 


for Example 5. 


Table 10: Comparison of absolute errors with other me 


10, 2 


2,71 


hods for Example 5 


tj present method Bernoulli Runge-Kutta Variational iteration 
method in [18] | method in [18] method in [18] 

0.2 0.0000 5.31 x 10-18 4.43 x 10-4 2.98 x 10-10 

0.4 | 1.0000 x 10724 1.24 x 10712 3.85 x 1073 1.01 x 1078 

0.6 0.0000 2.89 x 10~10 1.39 x 107? 8.24 x 10-8 

0.8 | 4.0000 x 10—2° 6.49 x 10710 3.53 x 1072 3.76 x 107-7 

1.0 | 1.0000 x 10-20 1.23 x 10-9 7.34 x 1072 1.26 x 107° 
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— exact solution © approximate solution — exact solution = y=1.6 =" y=L7 1 y=L8 —" y=19 


Figure 6: Comparison of (a) exact and approximate solutions, (b) u(t) for a = 8 = — 
and N=7,0= 7 and y = 1.6, 1.7, 1.8,1.9 for Example 6. 


i 
4? 


Example 6. Consider the following nonlinear fractional pantograph differ- 
ential equation of [18] : 


t 
D u(t) =1—- 2u*(5), 1<y7<2. 


The initial conditions are u(0) = 1, u’(0) = 0 and the exact solution is 
u(t) = cos(t) in case y = 2. By using the approximations introduced in the 
previous sections, we get the following algebraic equation: 


@OT (t) C4+2UT 0 OM (t)-1 20, U= Lise) pot C+LA ,y F, (35) 


where u(0) = 1 is approximated by ®(%)7(t) F and U is the operational ma- 
trix of product corresponding to the vector U. By choosing N = 7, equation 
(35) is collocated at roots of P{°’?)(«x). Then, by solving the resultant alge- 
braic system by means of Newton’s iteration method, the unknown vector C 
is determined for some values of a and (6 parameters. Table 11 displays the 
absolute errors for a = 8 = 3/2, y = 2, and various values of o. In Table 12, 
the values of the approximate solutions are shown at some arbitrary points 
for N = 7, a = 6B = —1/4, o = 1/2, and various values of 7. A graphical 
comparison between the exact and approximate solutions is found in part (a) 
of Figure 6 for a = 6 = —1/4, 0 = 1, and 7 = 2. The plots of u(t) have been 
compared in part (b) for y = 1.6,1.7,1.8,1.9, and o = 1/2. It is seen that 
the approximate solutions converge to the exact solution as y approaches 2. 
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Table 11: Absolute errors for various values of ¢, and N = 7 for Example 6 

t; o = 0.5 o=1 0 =1.5 o=2 

0.0 | 2.4053 x 10-7 | 6.0568 x 10-® | 5.6194 x 107° | 1.1748 x 107? 
0.2 | 4.0376 x 10-8 | 7.7581 x 1071! | 2.3175 x 107% | 1.0518 x 107? 
0.4 | 4.1630 x 1078 | 4.1069 x 107!° | 1.7587 x 1073 | 7.6767 x 1073 
0.6 | 8.1380 x 10-8 | 5.3919 x 10-1° | 9.8358 x 1074 | 3.3268 x 1073 
0.8 | 2.3864 x 10-7 | 2.3463 x 1071° | 3.8368 x 1074 | 2.4282 x 1073 
1.0 | 5.8922 x 10-® | 5.4779 x 10-9 | 4.2678 x 1074 | 9.5443 x 107% 


Table 12: Absolute errors for various values of y, and N = 7 for Example 6 


ti y= 1.6 y=1.7 y=1.8 = 1:9 

0.0 | 1.9622 x 10-7 | 1.1006 x 10-” | 2.1920 x 10-7 | 5.1996 x 1078 
0.2 | 1.1452 x 107® | 2.9183 x 1075 | 4.5715 x 1075 | 2.5122 x 1075 
0.4 | 1.4694 x 10-5 | 4.2206 x 10-5 | 7.3289 x 107° | 7.9829 x 1075 
0.6 | 1.4420 x 10-5 | 4.6407 x 10-5 | 8.8547 x 1075 | 1.0437 x 1075 
0.8 | 1.0683 x 10-5 | 4.2124 x 10-5 | 9.0686 x 107° | 1.1652 x 1074 
1.0 | 4.2390 x 10-® | 3.0183 x 10-® | 7.9795 x 10-® | 11529 x 1074 


9 Conclusion 


In this paper, the Jacobi operational method was used to numerically 
solve the linear and nonlinear fractional pantograph differential equations. 
This approach was based on the shifted fractional-order Jacobi functions 
on the interval [0,1], which converted a given linear or nonlinear fractional 
pantograph differential equation into a set of linear or nonlinear algebraic 
equations. Another advantage of the proposed method was that the coefhi- 
cients of the Jacobi series solution were easily found by using the presented 
algorithms. These algorithms decrease the computational costs compared 
to other methods such as the Bernoulli collocation method as it applies the 
lesser numbers of operational matrices. By using the introduced algorithms, 
the operational matrices of fractional integration, product, and pantograph 
were simply constructed. The illustrative examples confirmed the validity 
and simplicity of the proposed approach. It was found that the absolute er- 
ror of this method is lesser than the errors of other common methods and the 
errors of our proposed method were almost constant on the studied interval. 
This demonstrates the stability of the applied Jacobi collocation method. II- 
lustrative examples with satisfactory results were used to demonstrate the 
application of this method. Suggested approximations make this method 
very interesting and effective to the good agreement between approximate 
and exact solutions. Almost the results obtained by the technique developed 
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in this paper were more accurate than the results obtained by the Bernoulli 
collocation, Runge-Kutta, variational iteration, and 6- methods. 
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